Analyze the ABC output

Inputs

Plot formatting

Extract samples

Plot ABC posterior marginal distributions

$R_0$ posterior distributions

The $R_0$ value is given by $R_0 = \int^\infty_0 S_\gamma(t)\beta(t)dt$, where $S_\gamma(t)$ is the survival function of $\gamma$ and $\beta$ is the contact interval hazard. For a Weibull distribution, the survival function is $S_\gamma(t) = e^{-(t/\gamma_\theta)^{\gamma_\alpha}}$. Combining with the definition of $\beta$ we obtain
$R_0 = \int^\infty_0 e^{-(t/\gamma_\theta)^{\gamma_\alpha}}\frac{\beta_\alpha}{\beta_\theta}\left(\frac{t}{\beta_\theta}\right)^{\beta_\alpha-1}dt.$
We now estimate this integral by truncating it to a finite domain with error at most tol, and then use a trapezoidal approximation with npts many sample points.

Parameter Posterior CI's

Pearson correlation coefficients

Analyze Best Fit

Simulate the best fit

Plot equation coefficients

Plot best fit solution

Plot errors with ODH